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Abstract 

We implement for the first time the simulated annealing method (SAM) to the problem of 
searching for the global minimum in the hyper-surface of the chi-square function which depends on 
the values of the parameters of a Skyrme type effective nucleon-nucleon interaction. We undertake 
a realistic case of fitting the values of the Skyrme parameters to an extensive set of experimental 
data on the ground state properties of many nuclei ranging from normal to exotic ones. The set 
of experimental data used in our fitting procedure includes the radii for the valence W5/2 an d 
I/7/2 neutron orbits in the 17 and 41 Ca nuclei, respectively, and the breathing mode energies for 
several nuclei, in addition to the typically used data on binding energy, charge radii and spin-orbit 
splitting. We also include in the fit the critical density p cr and further constrain the values of the 
Skyrme parameters by requiring that (i) the quantity P = 3p4j|, directly related to the slope of 
the symmetry energy S, must be positive for densities up to 3po (ii) the enhancement factor k, 
associated with the isovector giant dipole resonance, should lie in the range of 0.1 — 0.5 and (hi) the 
Landau parameter G' is positive at p = pq. We provide simple but consistent schemes to account 
for the center of mass corrections to the binding energy and charge radii. 

PACS numbers: 21.10.Dr, 21.10.Pc, 21.30.Fe,21.60.Ka 



I. INTRODUCTION 

Since the pioneering work of Brink and Vautherin , continuous efforts have been made 
to readjust the parameters of the Skyrme type effective nucleon-nucleon interaction to better 
reproduce experimental data. Most of the parameters of the Skyrme interactions available 
in the literature were obtained by fitting the Hartree-Fock (HF) results to experimental data 
on bulk properties of a few stable closed shell nuclei. Only recently, several families of the 
Skym e p_ S e,, SkIl-6, SL yl - 10 , SKX m d Sk0 SflflflflQ were obtained b y 
fitting HF results to the experimental data on the bulk properties of nuclei ranging from the 
(3— stable nuclei to those near the proton and/or neutron drip lines. In the SKX interaction, 
to yield appropriately the values for the binding energy differences between mirror nuclei also 
referred to as the Coulomb displacement energy (CDE), the contribution of the Coulomb 
exchange term in the HF equations is ignored and the direct Coulomb term is evaluated 
by replacing the pointproton distribution by its charge distribution. However, it has been 



further shown hi Ref. Q thai the SKX iutomotiou is not suitable for studying tin* properties 
of neutron stars. Since, for the SKX interaction, the quantity 

dS , . 

p = 3 »t p > « 

which is directly related to the slope of the symmetry energy coefficient S, becomes negative 
for nuclear matter densities p well below 3po (Po — 0-16 fm~ 3 is the saturation density). On 
the other hand, the SkIf-6, SLyl-10 and SkO Skyrme interactions are found to be suitable 
for the study of neutron stars |8|. But, these families of Skyrme interactions significantly 
underestimate the values of the CDE for mirror nuclei, since the Coulomb exchange term 
was included. Thus, it is desirable to have a unified interaction which includes the merits of 
several families of the Skyrme interactions as mentioned above. One can further enhance the 
applicability of the Skyrme type effective nucleon-nucleon interaction by imposing certain 
constraints as discussed below. 

The aim of this work is twofold. We implement, for the first time, the simulated annealing 
method (SAM) ^ Q to fit the values of the Skyrme parameters and develop a more 
realistic Skyrme type effective interaction. The SAM is an elegant technique for optimization 
problems of large scale, in particular, where a desired global extremum is hidden among many 
local extrema. This method has been found to be an extremely useful tool for a wide variety 
of minimization problems of large non-linear systems in many different areas of science (e.g., 



see Refs. jll[ 12, 13]). Very recently 14, Q|, the SAM was used to generate some initial 
trial parameter sets for the point coupling variant of the relativistic mean field model. In 
the present context, we use the SAM to determine the values of the Skyrme parameters by 
searching for the global minimum in the hyper-surface of the x 2 function given as, 

X N d -N p ^\ a t J U 

where, and N p are the number of experimental data points and the number of fitted 
parameters, respectively, <Tj is the theoretical error and M^ xp and Mf 1 are the experimental 
and the corresponding theoretical values, respectively, for a given observable. The values of 
X 2 depends on the Skyrme parameters, since, the Mf 1 in Eq. (JSJ) is calculated using the HF 
approach with a Skyrme type effective nucleon-nucleon interaction. 

Toward the purpose of obtaining a more realistic parameterization of the Skyrme interac- 
tion, we apply the SAM to fit the HF results to an extensive set of experimental data for the 
binding energy, charge radii, spin-orbit splitting and root mean square (rms) radii of valence 
neutron orbits. Our data set used in the fit consists of 14 spherical nuclei, namely, 16 0, 24 0, 
34 Si, 40 Ca, 48 Ca, 48 Ni, 56 Ni, 68 Ni, 78 Ni, 88 Sr, 90 Zr, 100 Sn, 132 Sn and 208 Pb. We also include in 
our fit the experimental data for the breathing mode energies for the 90 Zr, 116 Sn, 144 Sm and 
208 Pb nuclei. In addition, we include in the fit, the critical density p cr which is determined 
from the stability conditions for the Landau parameters. r], G| • We further constrain the 
values of the Skyrme parameters by requiring that (i) the quantity P must be positive for 
densities up to 3po! a condition imposed by neutron star models [8j], (ii) the enhancement 
factor k, associated with the Thomas-Reiche-Kuhn (TRK) sum rule for the isovector giant 
dip* resonance, lies in the r a„ g e of 0, - 0.5 SEfflQ and (id) the Landau pa— 
G' , crucial for the spin properties of finite nuclei and nuclear matter, should be positive at 



P = Po 
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2l| . We also provide simple but consistent schemes to account appropriately for 
the CDE and the center of mass corrections to the binding energy and charge radii. In order 
to check the reliability of the proposed Skyrme interactions for the study of high density 
matter, we have examined in detail the behavior of the symmetry energy and the nature of 
the equation of state (EOS) for pure neutron matter at densities relevant for the neutron 
star models. 

We have organized our paper as follows. In Sec. |H] we briefly outline the form of the 
Skyrme nucleon-nucleon (NN) effective interaction and the corresponding energy density 



functional adopted in the present work. In this section, we also provide a viable strategies 
for the calculations of CDE and the center of mass corrections to the total binding energy 
and charge radii. In Sec. II 111 we provide the relations between the Skyrme parameters and 
the various nuclear matter properties, which we have used to implement the SAM algorithm 
as described in Sec. II VI The set of the experimental data along with the theoretical errors 
and the constraints used in the fit to determine the values of the Skyrme parameters are 
given in Sec. In Sec. IVII we present our results for the two different fits carried out in 
this work. Finally, in Sec. I VI II we summarize our main results and discuss the scope for the 
further improvement of the present work. 



II. SKYRME ENERGY DENSITY FUNCTIONAL 



In this work^we adopt the following form for the Skyrme type effective nucleon-nucleon 
interaction 
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where, t iy x,, a and Wo are the parameters of the interaction and P[ 2 is the spin exchange 
operator, ~~a i is the Pauli spin operator, k \ 2 — — z(Vi — V 2 )/2 and k \ 2 — — i( Vi — V 2 )/2 . 
Here, the right and left arrows indicate that the momentum operators act on the right and 
on the left, respectively. The corresponding mean-field Vhf and the total energy E of the 
system are given by, 



E = / H(r)d 3 r 



_5H 

where, the Skyrme energy density functional H(r), obtained using Eq. (jSJ), is given by 



(4) 

US, 



H = JC + Hq + Hr + H, 



cff 



H 



fin 



H so + H sg + Tt 



t-Coul 



(5) 



where, /C = -|-r is the kinetic energy term, Hq is the zero-range term, TC S the density 



dependent term, 7i eff an effective-mass term, 7i fin a finite-range term, TC so a spin-orbit term, 



H sg a term due to tensor coupling with spin and gradient and Hcoui is the contribution to 
the energy density due to the Coulomb interaction. For the Skyrme interaction of Eq. @, 
we have, 



H = ~i [(2 + x )p 2 - (2x + l)(pj + pi)) , (6) 



^3 = ^3P° [(2 + x 3 )p 2 - (2x 3 + 1) (P 2 P + Pi)} , (7) 



1 1 

H eS = - [h(2 + xt) + t 2 (2 + x 2 )} rp + - [t 2 (2x 2 + 1) - ti(2xi + 1)] (r pPp + r n p n ), (8) 

o o 

ftfin = — [3ti(2 + a;i) - t 2 (2 + x 2 )] (Vp) 2 

-1 [3*1(2*! + 1) + t 2 (2x 2 + 1)] [(Vp p ) 2 + (Vp n ) 2 ] , (9) 

Wn 

n so = -^[j- vp + x w (j p - vp p + j n - Vp n )], (io) 

Hsg = -ig(*lZl + t 2^ 2 )J 2 + — (tl - t 2 ) [Jp 2 + Jn'] • (11) 

Here, p = p p +p n , T = Tp+Tm and J = J p +J„ are the particle number density, kinetic energy 
density and spin density with p and n denoting the protons and neutrons, respectively. We 
have used the value of h 2 /2m = 20.734 MeVfm 2 in our calculations. We would like to 
emphasize that we have included the contributions from the spin-density term as given by 
Eq. (fTTjl which is ignored in many Skyrme HF calculations. Although the contributions 
from the Eq. (fTTf to the binding energy and charge radii are not very significant, they are 



very crucial for the calculation of the Landau parameter G' Q 



221. 



A. Coulomb energy 

The contribution to the energy density (Eq. (|SJ)) from the Coulomb interaction can be 
written as a sum of a direct and an exchange terms, 

Wo Ml (r) = 7^(0+7^(0. (12) 

For the direct term it is common to adopt the expression 

W&i(r) = \e 2 Pp {r) j (13) 



ion 0, 



and for the corresponding exchange term to use the Slater approximat 

^ oul (r) = -%(r)f^V /3 . (14) 
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It is very important to emphasize that the definitions of Eqs. f|13|) and ()14|) are not for 
the bona fide direct and exchange terms, since each of them includes the contributions of 
the self- interaction, which appear in opposite signs and cancel out in Eq. (|12|). Note, in 
particular, that the direct term of Eq. (J13|) is proportional to Z 2 and not to ZiZ — 1), as 
it should be for a direct term, see a detailed discussion in Ref. [24]. We point out that for 
the CDE of mirror nuclei the magnitude of the self interaction term is CDE/(2Z), i.e., one 
has a spurious increase in the calculated CDE of about 6.3% and 2.5% for the A=17 and 41 
systems of mirror nuclei, respectively. 

We recall that within the mean-field approximation, adjusted to reproduce the experi- 
mental values of the charge rms radii, the calculated CDE of analog states (obtained using 
Eq. ()12|) ) are smaller than the corresponding experimental values by about 7%. It was 
irst shown in Ref. 25j that this discrepancy, also known as the Nolen-Schiffer anomaly 
2fi| . can be explained by taking into account the contributions due to long range correla- 
tions (LRC) and due to the charge symmetry breaking (CSB) in the NN interaction, see 
also Ref. [221 ■ We add that for the mirror nuclei with A=17 and A=41, the LRC and the 
CSB each accounts for about half of the discrepancy between theory and experiment. Also, 
the magnitude of the bona fide exchange Coulomb term is about the same as that due to 
LRC, but with opposite sign. Therefore, neglecting the bona fide Coulomb exchange term 
does not resolved the discrepancy between theory and experiment, but can account for the 
contribution of LRC. We add that in Ref. it was shown that by ignoring the Coulomb 
exchange term in the form of of Eq. (14) in Eq. (12), i.e., by including only the Coulomb 
direct term in the form of Eq. (13) (as is the case for the SKX interaction), one reproduces 
the experimental values of the CDE. It should be clear that this is due to the fact that by 
adopting the form of Eq. (13) for the Coulomb direct term one not only neglects the bona 
fide Coulomb exchange term, but also adds the spurious contribution of the self-interaction 
term. The unphysical neglect of the bona fide Coulomb exchange term together with the 
spurious contribution of the self-interaction term results in a contribution to CDE which is 
similar in magnitude to that obtained from the LRC + CSB terms. For simplicity, we will 
adopt in this work the form of Eq. (13) for the Coulomb direct term. 



B. Center of mass corrections to the binding energy and charge radii 



The HF approach applied to finite nuclei violates the translational invariance, introducing 
a spurious center of mass (CM) motion. Thus, one must extract the contributions of the 
center of mass motion to the binding energy B, radii and other observables. An accurate way 
to restore the translational invariance is to use the projection method. But, it is numerically 
very expensive. So, it is desirable to develop simple schemes for the CM corrections to various 
observables. Normally, one makes the CM corrections only to the binding energy and not 
to the radii. However, the CM corrections to the rms radii for light nuclei may be as large 



as 2% 



2J]. In the present work we shall consider the CM corrections to the binding energy 



as well as charge rms radii used to fit the Skyrme parameters. 

To account for the CM correction to the total binding energy, one must subtract from it 
the so-called CM energy given as, 

Ecu = ^(P 2 ) (15) 

where, P = — ih X^li V* is the total linear momentum operator. Traditionally, one simplifies 
the computation of Eq. (j!5|) by taking into account only the one-body parts of it, which can 
be easily achieved by replacing — — > — [l — 4] , in the kinetic energy term. In this case, 
the effects of neglecting the two-body part of Eq. (|15j) are compensated by renormalization 
of the force parameters. This may induce in the forces an incorrect trend with respect to A 
which becomes visible in the nuclear matter properties. In fact, it is found in Ref. that 
an over simplified treatment of Eqm obtained by renormalizing the nucleon mass appearing 
in the kinetic energy term leads to a larger value of the surface energy coefficient than those 
obtained using the full CM correction. This gives rise to differences in the deformation 
energy which becomes quite pronounced for the super deformed states. Very recently jl^, 
we also find that a large value of the surface energy coefficient yield a smaller value for the 
critical density. Thus, an appropriate and still simple scheme to evaluate Eq. (|T5)l is highly 
desirable. We note, however, that the SLy6, SLy7 and SLylO interactions |4( have been 
obtained by evaluating Eq. (JTSJl ( i.e., including the one and two-body CM terms of Eq. 
(j!5Jl ). In the harmonic oscillator (HO) approximation, Ecm of Eq. (|15|) is given by, 

E°cm = (16) 
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A value of fkv — 41A" 1 / 3 MeV is used in many relativistic mean- field calculations 30 . 13I | . 
An improved version for the CM correction can be obtained by modifying the oscillator 
frequency as hw = — 25A~ 2 ^ 3 MeV, which has been used in Ref. |5[ to obtain the 

SKX interaction. Here, we employ a simple but more consistent scheme to evaluate the 
Ecm using the HO approximation. We determine the oscillator frequency hu appearing in 
Eq. (|16|) using the mean square mass radii (r 2 ) calculated in the HF approach as, 

tuv = ^— y* ^ + -], (17) 

mA(r 2 } ^ l 2 J K J 

where, the sum runs over all the occupied single-particle states for the protons and neutrons 
and Mi is the oscillator quantum number. We emphasize that this scheme is quite reliable 
even for the nuclei away from the /3-stable line where the values of the rms radii deviate from 
the A 1 / 3 law. We have calculated the total binding energy for the SLy7 interaction using our 
simple scheme for the CM correction, Eq. (|17|h and compare them with those given in Ref. 
4]], obtained by using one-body and two-body parts of the Eq. (|15|) . For example, we find 
that for the 16 0, 40 Ca, 132 Sn and 208 Pb nuclei the total binding energy B = 128.65 (128.55), 
344.98 (344.90), 1102.38 (1102.77) and 1636.29 (1636.76) MeV, respectively, where the values 
in parenthesis are taken from Ref. |4|. This clearly indicates that the CM correction to the 
binding energy can be reliably estimated using Eq. ()17j) . We would also like to remark, 
however, that the Ecm calculated using the oscillator frequency as frw = 45A _1//3 — 25A~ 2 / 3 
MeV in Eqs. (J 16)) overestimates the value of binding energy in light nuclei (e.g., 16 and 
40 Ca) by about 1 — 2 MeV which is quite significant. 

The mean-square radius for the point proton distribution corrected for the CM motion 
is obtained as |24j |. 

(4) = (rJ)HF - (18) 

where, v = muo/h is the size parameter. Therefore, the corresponding mean-square charge 
radius to be fitted to the experimental data is obtained as, 



(r 2 ch ) = (r 2 p )n F - ^4 + ( r2 ) P + f (^)n + z 



\{— )£cy +1 )M>- z >y ( 19 ) 

/> \mc ) z — > 



nljr 

where, (r 2 ) p and (r 2 ) n are the mean-squared radii of the proton and neutron charge distri- 
butions, respectively. The last term in Eq. (fT9")l is due to the spin-orbit effect j^. We use, 
(r 2 ) n = -0.12 fm 2 and the recent 3 value of ( r2 ) P = °- 801 im<2 - 
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C. Determination of the critical density 



We use the stability conditions of the Landau parameters for the symmetric nuclear 
matter and pure neutron matter to calculate the critical density p cr for the Skyrme type 
effective nucleon-nucleon interactions. The stability conditions are given as 3^1 , 

^ > -{21 + 1), (20) 

where, Ai stands for the Landau parameters Fi, F(, Gi and G\ for a given multipolarity I. 
Skyrme interactions only contain monopolar and dipolar contributions to the particle-hole 
interaction so that all Landau parameters are zero for I > 1. Thus, there are 12 different 
Landau parameters, i.e., Fi, F{, Gi and G\ (I = 0, 1) for the symmetric nuclear matter and 
G\ n) (I = 0,1) for the pure neutron matter. Each of these parameters must satisfy 
the inequality condition given by Eq. ()20)1 . Explicit expressions for the Landau parameters 
in terms of the Skyrme parameters can be found in Refs. Q, 0|. The critical density 
is nothing but the maximum density beyond which at least one of the Landau Parameter 
does not satisfy Eq. (|20p. Following Ref. r|, one can obtain the values of the Landau 
parameters at any density for a given set of the Skyrme parameters. Thus, for a given set 
of Skyrme parameters one can easily obtain the value of p cr . As mentioned in Sec. HJ we 
include p cr in the fit. 

D. Breathing mode energy 

We also include in our fit the experimental data on the breathing mode energy for several 
nuclei. We consider the fully self-consistent values for the breathing mode constrained energy 
defined as, 

* = J^, (21) 
V m_i 

where are the energy moments 

/■oo 

m k = u k S{u)du, (22) 
Jo 

of the strength function 

5M = ^|(n|F|0)| 2 5( W -^ n ), (23) 



for the monopole operator F(r) = ^;=i f( r i), with /(r) = r 2 . The moments m^ for k — — 1 
and 1 appearing in Eq. (|2*T|) can be obtained using the constrained HF (CHF) and the 



double commutator sum rule, respectively 
in terms of the ground state density p as, 



37 



3a ]. The moment mi can be expressed 



mi 



fi 2 

2-(r 2 ), 
m 



(24) 



where, 

As described in detail in Ref. 
is given as, 



36 
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(r 2 ) = / r 2 p(r)dr. 



(25) 

381 ] . m_i can be evaluated via the CHF approach and 
1 d 



2d\ 



(rl) 



(26) 



A=0 



where, (r 2 ) = ($a \r 2 \ $a), with <3>a being the HF solution to the constrained HF Hamiltonian 
H-Xf. 

III. SKYRME PARAMETERS AND NUCLEAR MATTER PROPERTIES 

In this section we discuss the relationship between the Skyrme parameters and the various 
quantities describing the nuclear matter. In the next section we use these relations to 
implement the SAM algorithm. The Skyrme parameters ti, Xi and a for a fixed value of Wq 
can be expressed Jin terms of the quantities associated with the symmetric nuclear matter 
as follows 3], 16 



(-B/A + (2m/m* - 3) (fi 2 /10m) kf) (±K nm - (1 - 6m* /5m) (h 2 /9m*) kj^ 
-B/A + \K nm - (4m/3m* - 1) (h 2 /10m) k 2 ^ 



(4m/3m* - 1) (^ 2 /10m) k 2 



+ 1- 
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h = - [T + T s 



2m* k 
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K 2 TC 2 
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m Pnm 
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p^t 1 -B/A + \K nm - (4m/3m* - 1) (H 2 /Wm) kf 



(28) 
(29) 
(30) 
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(32) 



x 2 = — [8T - 3ta - 5t 2 ] 

4t2 



(33) 



where, 



Ct+l 
nm 



a 



' h 2 l 

— fc? - — ((4 + 5x 2 )t 2 - 3tiXi) p nm k 2 f -3J + L 
bm J 12 1 

B/A - \K nm + (4m/3m* - 1) (h 2 /10m) k 2 
-B/A + (2m /m* - 3) (ft 2 /10m) fc2 ' 



(34) 



(35) 



^0=o ( 3t i + ( 5 + 4x M 



H 2 



mp ri 



m 



rrr 



- 1 



(36) 



and 



T s = \ [9*i - (5 + 4a; 2 )t 2 ] 



3^^ 



1/3 



(37) 



(38) 



In Eqs. (J27J) - (|35p. the various quantities characterizing the nuclear matter are the binding 
energy per nucleon B/A, isoscalar effective mass m*/m, nuclear matter incompressibility 
coefficient K nm , symmetry energy coefficient J = S(p = p nm ), the coefficient L = P(p = 
Pnrn)-, enhancement factor k and Landau parameter G' . All these quantities are taken at 
the saturation density p nm . It must be pointed out that the expression for the parameter 
G'q used in Eq. (J29|) includes the contributions from the spin-density term present in the 
Skyrme energy density functional |22j]. So, for consistency, the HF calculations are also 
performed by including the contributions from the spin density. Once, To is known, T s can 
be calculated for a given value of the surface energy E s as jl^ . 



E s = 87tr z / dp 
'o 



h 2 5 1 m* 2 

3^-3^+8^-7^ 



1/2 



[B(p nm )/A-B(p)/A] 1/2 , (39) 
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where, B(p)/A is the binding energy per nucleon given by, 



B(p) 
A 



and, 



r 



1 1/3 



4:%p nn 



9 
16 



Wo 2 . 



(40) 

(41) 
(42) 



The manner in which Eqs. (127(1 - (|33|l can be used to evaluate the Skyrme parameters tj, 
Xi and a is as follows. At first, the parameters to and a can be calculated in terms of B/A, 
p nm , K nm and m*/m, using Eqs. ((77(1 and ((33(1 . Then, the parameter t 3 can be determined 
using Eq. (|5U|l. Next, T and T s can be calculated using Eqs. ((Ho]) and (|39j). respectively. 
Once, the combinations To and T s of the Skyrme parameters are known, one can calculate 
the remaining parameters in the following sequence, ti, £2, £2, £1, ^3 and xq. 



IV. SIMULATED ANNEALING BASED ALGORITHM FOR THE MINIMIZA- 
TION OF x 2 

The simulated annealing method (SAM) is a generalization of a Monte Carlo technique, 
based on the Metropolis algorithm |40(, initially developed for examining the equation of 
the state of a many body system. The concept of SAM is based on the manner in which 
liquids freeze or metals recrystallize in the process of annealing. In an annealing process a 
metal, initially at high temperature and disordered, slowly cools so that the system at any 
time is in a thermodynamic equilibrium. As cooling proceeds, the system becomes more 
ordered and approaches a frozen ground state at zero temperature. 

With this brief background, we shall now implement the SAM to search for the global 
minimum of x 2 function as given by Eq. (j2J). One of the crucial key ingredients required 
to implement the SAM, in the present case, is to specify the lower and the upper limits for 
each of the Skyrme parameters. So that the global minimum for the x 2 is searched within 
these limits. However, from the literature (e.g., see Refs. BQ) 

we find that the Skyrme 

parameters vary over a wide range. To make the search process more efficient, we make use 
of the fact that most of the Skyrme parameters can be expressed in terms of the various 
quantities related to the nuclear matter properties as described in Sec. IIHI Most of these 
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nuclear matter quantities are known empirically within 10% — 20%. For convenience, we 
define a vector v with 10 components as, 

v = (B/A, K nm */m,E„J,L,K,G' ,W ). (43) 

Once the vector v is known we can calculate the values of all the Skyrme parameters as 
discussed in Sec. IIHI We also define the vectors v , Vi and d. The vector v and Vi contains 
the lower and the upper limits of each of the components of the vector v. The vector d 
represents the maximum displacement allowed in a single step for the components of the 
vector v. We implement the SAM algorithm using the following basic steps, 

(i) We start with a guess value for the vector v and calculate \ 2 ( sa Y) Xoid) usrn g Eq. (J2J) 
for a given set of the experimental data and the corresponding HF results together 
with the theoretical errors. 

(ii) We generate randomly a new set of Skyrme parameters using the following steps. First, 
we use a uniform random number to select a component v r of the vector v. Second, 
the randomly selected component v r is then assigned a new value, 

V r — > V r + T]d r , (44) 

where rj is a uniform random number which lies within the range of —1 to +1. The 
second step is repeated until the new value of v r is found within its allowed limits 
defined by vo and vi. We use this modified v to generate a new set of Skyrme 
parameters. It may be noted that a change in the value of a component of the vector 
v may lead to changes in the values of several Skyrme parameters. For example, a 
change in the value of K nm will alter the values of the Skyrme parameters to, £3 and 
a. 

(iii) The newly generated set of the Skyrme parameters is accepted by using the Metropolis 
algorithm as follows. We calculate the quantity, 

V{ X 2 ) = e (x °id-xLw)/T ; ( 45 ) 

where xLm * s obtained by using the newly generated set of the Skyrme parameters 
and T is a control parameter (an effective temperature). The new set of Skyrme 
parameters is accepted only if, 

V( X 2 ) > P, (46) 
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where (3 is a uniform random number which lies between and 1. If the new Skyrme 
parameters are accepted (i.e. Eq. (jjfjj) is satisfied), it is called a "successful reconfig- 
uration" . 

To search for the global minimum of \ 2 we begin with some reasonable value of an 
effective temperature T = Tj. For a given T iy we repeat steps (ii) and (iii) for, say, 100 N p 
reconfigurations, or for 10N P successful reconfigurations, whichever comes first. Then, we 
reduce the temperature by following a suitable annealing schedule. One encounters various 
annealing schedule available in the literature such as linear, exponential, Boltzmann and 
Cauchy annealing schedules Among these, the Boltzmann annealing schedule is the 

slowest one and the exponential annealing schedule is the fastest one. In the present work 
we have employed the Cauchy annealing schedule given by, 

T(k) = T t /ck (47) 

where, c is a constant, which is taken to be unity in the present work, and k — 1,2, 3, is 

the time index. We keep on reducing the value of T using Eq. (|37|) in the subsequent steps 
until the effort to reduce the value of x 2 further becomes sufficiently discouraging. 

In Table H] we list the values of all the components of the vectors v, vo, Vi and d used 
in the numerical computation. We have varied the components of the vector v over a wide 
range. The values of the maximum displacement as defined by the components of d are so 
chosen that the corresponding component of the vector v can be varied over the entire range 
given by the vectors vo and vi, within the adopted number of reconfigurations. We have 
carried out several sample runs and found that Tj = 1.25 along with the Cauchy annealing 
schedule yields reasonable values of the Skyrme parameters. We must mention here that the 
range for the quantities L, k and G' as given in Table H] are so chosen that they vary within 



acceptable limits 



V. EXPERIMENTAL DATA AND SOME CONSTRAINTS 

In this section we discuss our selection of the experimental data and the corresponding 
theoretical errors adopted in the \ 2 fit, Eq. (J2J), to the HF results. In Table ITT1 we summarize 
our choice of the experimental data. It must be noted that in addition to the typically used 
data on the binding energy, charge radii and spin-orbit splitting, we also include in our fit the 
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experimental data for the radii of valence neutron orbits and the breathing mode energies of 



42 



43, 



44, 



45 



46 



several nuclei. All of these experimental data are taken from Refs 
For the binding energy we use in our fit the error of 1.0 MeV except for the 100 Sn nuclei. 
The binding energy for the 100 Sn nucleus is determined from systematics and are expected 
to have large errors. Thus, we assign them a theoretical error of 2.0 MeV. For the charge 
rms radii we use the theoretical error of 0.02 fm except for the case of 56 Ni nucleus. The 
charge rms radii for the 56 Ni nucleus is obtained from systematics and we use the theoretical 
error of 0.04 fm. We consider in our fit the experimental data for the spin-orbit splittings 
for the 2p neutrons and protons in the 56 Ni nucleus and the rms radii for the 1 0^5/2 and I/7/2 
neutron orbits in 17 and 41 Ca nuclei, respectively. We use 

{1.88 MeV Neutrons 
(48) 
1.83 MeV Protons 

where, e is the "bare" single-particle energy obtained by unfolding the experimental data 
for the energy levels in 57 Ni and 57 Cu nuclei by appropriately accounting for the coupling 
to excitations of the core. Of course, it is more appropriate to use the splitting of high I 
orbits in a heavy nucleus (e.g., 208 Pb nucleus) to determine the strength of the spin-orbit 
interaction. But, to the best of our knowledge, unlike for the 56 Ni nucleus the bare single- 
particle energies for the heavier nuclei are not available. For the rms radii of the valence 
neutron orbits in 17 and 41 Ca nuclei we use r„(z/lc?5/ 2 ) = 3.36 fm and r„(z/l/ 7 / 2 ) = 3.99 fm, 



44. 



45j respectively. The theoretical error taken for the spin-orbit splitting data is 0.2 MeV 
and for the rms radii for the valence neutron orbits we use the experimental error of 0.06 fm. 
We must point out that the choice of the theoretical error on the rms radii for the valence 
neutron orbits is due to the large uncertainties associated with their extraction from the 
experimental measurements. To be consistent with the way these valence neutron radii are 
determined, we do not include the center of mass correction to these data. The experimental 
data for the breathing mode constrained energies E Q included in our fit are 17.81, 15.90, 



15.25 and 14.18 MeV for the 90 Zr, 116 Sn, 144 Sm and 208 Pb nuclei |42|, respectively, with the 
theoretical error taken to be 0.5 MeV for the 90 Zr nucleus and 0.3 MeV for the other nuclei. 
We also include the critical density p cr in the fit assuming a value of 2.5po with an error of 
0.5po- Further the values of the Skyrme parameters are constrained by requiring that (i) 
P ^ for p ^ 3p , (h) k = 0.1 - 0.5 and (iii) G' ^ at p = p . 
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VI. RESULTS AND DISCUSSIONS 



In the preceding sections we have described in detail the implementation of the SAM based 
algorithm to fit the values of the Skyrme parameters to a given set of the experimental data 
considered in this work. We have carried out two different fits. These fits are carried out by 
using the same set of experimental data along with some constraints as discussed in Sec. El 
We name these fits as, 

1. KDEO; only the direct Coulomb term in the form of Eq. (JT3J) is included. 

2. KDE; the direct as well as the exchange Coulomb terms are included (Eqs. ()12 )) -()14 p ). 

The CM corrections to the total binding energy, Eqs. (fTfij) and (jl7|) . and the charge rms 
radii, Eqs. (JTHJ) and (fTDJ) . are carried out using the schemes described in Sec. Ill Bl 

We shall first consider some technical aspects required to implement the SAM. As it is 
evident from Sec. IIVI there are two crucial ingredients, namely, (i) initial value for the control 
parameter T = Tj and (ii) annealing schedule which determines the subsequent value for T. 
These ingredients essentially controls the computer time and the quality of the final fit. If 
one starts with a smaller value for T« and/or uses a faster annealing schedule, one may not 
be able to hit the global minimum of the objective function and rather get trapped in one of 
the local minima. In the present work we have employed the Cauchy annealing schedule. We 
have carried out several trial calculations and find that Tj = 1.25, along with the Cauchy's 
annealing schedule as given by Eq. (07J), yields reasonable values for the best fit parameters. 
To validate the present approach we carried out the following checks. Starting with the 
final values of the Skyrme parameters obtained using the SAM, we attempted to minimize 
further the value of \ 2 using the Levenberg-Marquardt (LM) method Q| as conventionally 
used. But, we found no further decrease in the value of the \ 2 - As an illustration, we plot 
in Fig. ^the average value (x 2 )t as an inverse function of the control parameter T for the 
KDEO case. The curves labeled v and v\ represent the results obtained from two different 
choices of the starting values for the Skyrme parameters. The initial value of the Skyrme 
parameters for the curve labeled v (solid line) and v\ (dashed line) are obtained using the 
set of values given in the second and fourth columns of Table respectively. The value of 
(x 2 )t is obtained by averaging over all the successful reconfigurations for a given T. We see 
from Fig. QJthat the value of (x 2 )t show a remarkable decrease at initial stages and then 
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oscillates before saturating to a minimum value for T ^ 0.05. The value of x 2 & t lower T is 
more or less independent of the starting values for the Skyrme parameters. In Fig|2]we have 
displayed the variation of Ax^ — ({x 2 ~ (x 2 )) 2 )t as an inverse function of T. We see that the 
fluctuations in the value of \ 2 is large for larger values of T. As T decreases, fluctuations in 
the value of x 2 a l so decrease rapidly. This means that initial value for T should not be too 
small, because, at smaller T it is less likely to jump from a configuration with lower value of 
X 2 to one having higher value. As a result, one may get trapped in a local minima. In Table 
IIHI we give the values of the parameters for the KDEO interaction at the minimum value of 
the x 2 obtained from different choices for the starting values for the Skyrme parameters. It 
is interesting to note that not only the final value of the x 2 is l ess sensitive to the choice 
of the initial parameters, but, also the resulting Skyrme parameters are also quite close. In 
what follows, we shall present the results for the KDEO and KDE interactions. The starting 
(or guess) values for he Skyrme parameters used to generate these interactions are obtained 
from the nuclear matter quantities given in the second column of the Table H] 

In Table HVl we give the values for the various quantities characterizing the nuclear matter 
obtained at the minimum value of the x 2 ■ We also note that the values of all the nuclear 
matter properties for the KDEO and KDE Skyrme interaction are closer to those obtained 
for the SLy7 interaction. However, it is worth mentioning that the values of the K nm and 
m*/m for both the interactions generated here emerge from the fit, unlike the SLy type 
interactions where the values for these quantities were kept fixed. In our fits, the values of 
the K nm and m*/m are mainly constrained by the inclusion of the experimental data on 
breathing mode energy and the value of critical density p cr = 2.5po ± 0.5po Q, l3- I n the 
last row of this table we give the values of x 2 & t the minimum. For the sake of completeness, 
we list in Table IVl the values of the Skyrme parameters obtained in the fits. One can easily 
calculate the values of these Skyrme parameters using the various nuclear matter quantities 
given in Table HVl as described in Sec. IHII In Table fVl we also give in parenthesis the values 
of the standard deviations for the Skyrme parameters. Since, within the SAM algorithm 
one can not calculate these standard deviations in a straight forward manner, we resort 
to some alternative approach. We have determined the values of the standard deviations 
on the parameters for the KDEO and KDE interactions using the LM method. The LM 
method requires two inputs, namely, set of the experimental data and the starting values of 
the interaction parameters. The set of experimental data is taken to be exactly the same 
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as the one used to generate the KDEO and KDE interactions. The starting values of the 
interactions parameters used are the ones obtained using SAM for the KDEO and KDE 
interactions. 

In Table IVT1 we present our results for the deviation AB = B exp — B th for the values of the 
binding energy obtained from the newly generated KDEO and KDE interactions. Similar de- 
viations Ar c h for the charge rms radii are presented in Table IVTT1 For comparison, in the last 
columns of these Tables we give the values of AB and Ar c h for the SLy7 interaction, taken 
from Ref. ^J. One can easily verify from Table IVD that the magnitude of the deviations for 
the binding energy for most of the cases is much less than 0.5% in case of KDEO interaction. 
The KDE interaction yields larger error in the values of the binding energy (~ 0.6 — 1.0%) 
for the 16 0, 48 Ni and 100 Sn nuclei. We would also like to remark here that in determining 
the SKX interaction, the binding energy for the 56 Ni nucleus was not considered in the fit 
and that for the 100 Sn nucleus was included in the fit with the theoretical error of 1.0 MeV. 
We find that if one attempt to do so, the binding energy for the 56 Ni becomes off by more 
than 3 MeV. We see from Table IVTll that, except for the 16 and 48 ca nuclei, the deviations 
in the values of the charge rms radii for the KDEO interaction is less than 0.5%. In addition 
to the binding energy and the charge rms radii of the nuclei used in our fits, we have also 
considered a few more experimental data as discussed in Sec. In Tables IVIIII and [D^ 
we present our results for these additional quantities. The values of p cr is greater than 2p - 
The values for the radii of valence neutron orbits and the spin-orbit splittings considered 
in our fits are quite reasonable for all the interactions considered here. It can be seen from 
Table IIXI that our fit to the breathing mode constrained energies are overall in reasonable 
agreement with the corresponding experimental data. 

We now consider our results for the binding energy difference between the 48 Ca and 48 Ni 
mirror nuclei. One may verify from Table [VTl that the binding energy difference S( 48 Ca) — 
5( 48 Ni) = 67.23 and 64.02 MeV for the KDEO and KDE interactions, respectively, compared 
to the experimental value of 68.85 MeV. We would also like to add that the said difference 
for the SKX interaction is 66.3 MeV which is about 1.0 MeV lower than our most realistic 
KDEO interaction. On the other hand, most of the Skyrme interactions which include the 
contribution from the exchange Coulomb term yield _B( 48 Ca) — 5( 48 Ni) pa 63 MeV, which 
is about 6 MeV lower than the corresponding experimental value. 
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We present in Table iXl our results for the neutron skin, r n — r p , the difference between the 
rms radii for the point neutrons and protons density distributions, for the KDEO and KDE 
interactions. We compare in Tables IXil and l_X.ll! the values of the single-particle energies 
with the available experimental data for the 40 Ca and 208 Pb nuclei |48l 149 1 . respectively. We 
find that the single-particle energies, for the occupied states near the fermi-energy compare 
reasonably well with the experimental ones. We would like to remark here that the HF 
approach alone is not expected to reproduce the experimental single-particle energies and 
there fore we have not included them in our fit. 

Finally, we consider the behavior of the symmetry energy coefficient S(p) for densities 
relevant to the study of neutron stars. It is well known |50l l51 | that the values of S(p) 
and the resulting EOS for pure neutron matter at higher densities (p > 2p ) are crucial in 
understanding the various properties of neutron star. For example, the proton fraction at 
any density depends strongly on the value of S(p) at that density, which in turn affects the 
chemical compositions as well as the cooling mechanism of the neutron star Yet, no 
consensus is reached for the density dependence of S(p). We display in Fig. El our results for 
the variation of the symmetry energy S as a function of the nuclear matter density p. We see 
for the KDEO and KDE interactions that the value of S increases with density for p < 3/?o- 
All of these interactions are quite suitable for modeling the neutron star with masses close 
to the canonical one Q], because, they yield S > for p < 4p . In Fig. HJwe plot the EOS 
for the pure neutron matter resulting from the KDEO and KDE interactions and com par e 
them with the ones obtained for SLy7 interaction and the realistic UV14+UVII model [53]. 
It is striking to note that our results for the KDEO and KDE interactions are in harmony 
with the EOS for the UV14+UVII model, though, unlike the SLy7 interaction we did not 
include in our fit the neutron matter EOS of the realistic UV14+UVII interaction. This 
seems to be due to the constrain imposed on the quantity P, which is related to the slope 
of the symmetry energy coefficient (see Eq. ) . 



VII. CONCLUSIONS 



We have implemented the simulated annealing method to fit the values of the parameters 
of the Skyrme interaction of Eq. (JHJ) by searching for the global minimum in the hyper surface 
of the x 2 function, Eq. To demonstrate the applicability of this method we have fitted 
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the values of the Skyrme parameters to an extensive set of experimental data together with 
a few additional constraints. Our experimental data set consists of the binding energies for 
14 nuclei ranging from the normal to exotic (proton or neutron rich) ones, charge rms radii 
for 7 nuclei, spin-orbit splittings for the 2p proton and neutron orbits of the 56 Ni nucleus and 
rms radii for lds/2 and I/7/2 valence neutron orbits in the 17 and 41 Ca nuclei, respectively. 
We also include in the fit the critical density p cr determined from the stability conditions 
for the Landau parameters. The additional constraints imposed on the Skyrme parameters 
are (i) the quantity P = 3p^, directly related to the slope of the symmetry energy S, must 
be positive for densities up to 3po; a condition imposed by the neutron star models p], (ii) 
the enhancement factor k, associated with the TRK sum rule for the isovector giant dipole 
resonance, should lie in the range of 0.1 — 0.5 and (iii) the Landau parameter G' , crucial for 
the spin properties of finite nuclei and nuclear matter, should be positive at p = po- 

Using these experimental data along with the additional constrains, we have carried out 
two different fits named as KDEO and KDE, as described in Sec. IVII The corrections to the 
binding energy and charge rms radii due to the center of mass motion were performed using 
simple but consistent schemes. The nuclear matter properties for both interactions proposed 
in the present work are obtained directly from the fit. The selection of the experimental 
data in conjugation with some constraints ensures that these interactions can be used to 
study the bulk ground state properties of nuclei ranging from the stable to the ones near the 
proton and neutron drip lines, as well as the properties of neutron stars. The interactions 
obtained in the present work encompasses the merits of the SKX and SLy type of Skyrme 
interactions. 

Before closing, we would like to mention that the method as well as the fitting strategy 
presented in this work can be improved in several ways. The SAM is a very adaptive 
approach and therefore it offers a significant scope for further improvement. For example, 
in the present work we jump from one configuration to another by randomly selecting a 
component of the vector v as defined by Eq. (|4*H|) . This selection was done using a uniform 
random number. However, one can think of performing random selection of a component 
of v by assigning a more plausible weight factors to these components. One can also try 
out different annealing schedules to determine the rate of cooling. In the present work we 
employed the Cauchy annealing schedule which yields a faster cooling rate than that of the 
Boltzmann schedule, but, a slower rate than the exponential annealing schedule. The effects 
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on the binding energy and radii due to the correlations beyond mean- field |5J, |55l |56| can 
be included in the fit. These effects are in particular important for the light nuclei. One 
may also include in the spin-orbit splitting the contributions due to the electromagnetic 
spin-orbit interacti on |46lj and modify the spin-orbit interaction by using the form proposed 
by Sagawa in Ref. [57fl. Last but not least, one may also include the experimental data on 
the giant dipole and quadrupole resonances while fitting the Skyrme parameters in addition 
to the breathing mode energy, as was done in the present work. 
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Figure captions 

FIG. 1: Variation of the average value of chi-square, (x 2 )t> as a function of the inverse of the control 
parameter T for the KDEO interaction for the two different choices of the starting parameters (see 
text for detail). 

FIG. 2: Variation of the fluctuations in the value of x 2 as a function of l/T for the KDEO 
interaction for the two different choices of the starting parameters (see text for detail). 

FIG. 3: Variation of the symmetry energy coefficient S(p) as a function of the nuclear matter 
density p. 

FIG. 4: Energy per particle for pure neutron matter as a function of density. Results for 

the two newly generated Skyrme interactions KDEO and KDE are compared with those obtained 
for the SLy7 Skyrme force and the realistic UV14+UVII model of Wiringa et al. 53]. 
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TABLE I: Values of the components of the vectors v, Vq, Vi and d used for implementing the 
SAM based algorithm for searching the global minimum of x 2 • The vector v initializes the value 
of x 2 , whereas, vo and vi limits the search space for the Skyrme parameters. The components of 
the vector d correspond to the maximum displacements allowed for the reconfiguration. 
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TABLE II: Selected experimental data for the binding energy B, charge rms radius r c h, rms radii 
of valence neutron orbits r v , spin-orbit splitting S-O, breathing mode constrained energy E$ and 
critical density p cr used in the fit to determine the parameters of the Skyrme interaction. 



Properties 


Nuclei 


Ref. 


B 


16 < 24 0, 34 Si, 40 < 48 Ca, 48,56,68,78^ 88g r> 90^ 100,132^ 208p b 


[41] 


r c h 


16 0, 40 - 48 Ca, 56 Ni, 88 Sr, 90 Zr, 208 Pb 


[42,43] 


r v (uld 5/2 ) 


17 


f44] 




41 Ca 


[45J 


S-0 


2p orbits in 56 Ni 


f46] 


E 


90 Zr, 116 Sn, 144 Sm, 208 Pb 


f47] 


Per 


nuclear matter 


see text 
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TABLE III: Comparison of the parameters for the KDEO interaction at the minimum value of \ 2 
obtained from different choices for the starting values of the Skyrme parameters. 



Parameter 


KDEO(v) 


KDEO(vi) 


f n fMeV-fm 3N l 

±V-1C< V 1111 i 


-2526 51 


-2553 08 


U fMeV-fm 5 ) 

\j 1 V IVllj V 1111 J 


430.94 


411.70 


fofMeV-fm 5 "! 

(/^ UvlC V 1111 } 


— 398 38 


—41 9 87 


t 3 (MeV-fm 3 ( 1+Q ) ) 


14235.5 


14603.6 




7583 


6483 


X\ 


-0.3087 


-0.3472 


X2 


-0.9495 


-0.9268 


X3 


1.1445 


0.9475 


^(MeV-fm 5 ) 


128.96 


124.41 


a 


0.1676 


0.1673 



TABLE IV: Nuclear matter properties for the KDEO and KDE interactions at the \ 2 = Xmi 



Parameter 


KDEO 


KDE 


SLy7 


B/A (MeV) 


16.11 


15.99 


15.92 


K nm (MeV) 


228.82 


223.89 


229.7 


Pnm 


0.161 


0.164 


0.158 


m* /m 


0.72 


0.76 


0.69 


E s (MeV) 


17.91 


17.98 


17.89 


J (MeV) 


33.00 


31.97 


31.99 


L (MeV) 


45.22 


41.43 


47.21 


K 


0.30 


0.16 


0.25 


G' 


0.05 


0.03 


0.04 




1.3 


2.2 
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TABLE V: The values of the Skyrme parameters for KDEO and KDE interactions obtained by 
minimizing the x 2 ■ For the sake of comparison we have also listed the values of the parameters for 
the SLy7 interaction. The values in Parenthesis are the standard deviation for the corresponding 
Skyrme parameters. 



Parameter 


KDEO 


KDE 


SLy7 


i (MeV-fm 3 ) 


-2526.51 (140.63) 


-2532.88 (115.32) 


-2482.41 


h (MeV-fm 5 ) 

J- V / 


430.94 (16.67) 


403.73 (27.63) 


457.97 


t 2 (MeV-fm 5 ) 


-398.38 (27.31) 


-394.56 (14.26) 


-419.85 


t 3 (MeV-fm 3 ( 1+Q )) 


14235.5 (680.73) 


14575.0 (641.99) 


13677.0 


x 


0.7583 (0.0655) 


0.7707 (0.0579) 


0.8460 


x\ 


-0.3087 (0.0165) 


-0.5229 (0.0298) 


-0.5110 


X2 


-0.9495 (0.0179) 


-0.8956 (0.0270) 


-1.0000 


X3 


1.1445 (0.0862) 


1.1716 (0.0767) 


1.3910 


W (MeV-fm 5 ) 


128.96 (3.33) 


128.06 (4.39) 


126.00 


a 


0.1676 (0.0163) 


0.1690 (0.0144) 


0.1667 
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TABLE VI: Results for the total binding energy B (in MeV) for several nuclei. The experimental 
data B exp used to fit the Skyrme parameters were taken from The theoretical error a was 
taken to be 2.0 MeV for the 100 Sn nucleus and 1.0 MeV for the other nuclei. In 3rd and 4th columns 



we give 


the values for AB = B exp 


— B th obtained from our new fi 


;s. The last column contains the 


values for AB for the SLy7 Skyrme interaction taken from Ref. ; 


I 








AB = B exp - B th 




in uciei 




KDEO 


KDE 


QT ir7 


I6n 


1 97 R90. 
1Z / .OZU 


0.394 


1.011 


— u.yo 


v > 


100.054 


-0.581 


0.370 




34c: 
ol 


Zoo.4fcZ / 


-0.656 


0.060 






o4Z. U0U 


0.005 


0.252 


9 

— Z.OO 




/lie; QQn 

410. yyu 


1 88 

\J . -L OO 


1 165 


U.ll 


48 Ni 


347.136 


-1.437 


-3.670 




56 Ni 


483.991 


1.091 


1.016 


1.71 


68 Ni 


590.408 


0.169 


0.539 


1.06 


78 Ni 


641.940 


-0.252 


0.763 




88 Sr 


768.468 


0.826 


1.132 




90 Zr 


783.892 


-0.127 


-0.200 




100 Sn 


824.800 


-3.664 


-4.928 


-4.83 


132 Sn 


1102.850 


-0.422 


-0.314 


0.08 


208p] 3 


1636.430 


0.945 


-0.338 


-0.33 
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TABLE VII: Results for the charge root mean square (rms) radii r c h (in fm). The experimental 
data used in the fit to determine the values of the Skyrme parameters are taken from Refs. [l] . 4^ |. 
The theoretical error a taken to be 0.04 fm for the 56 Ni nucleus and 0.02 fm for the other nuclei. 
In 3rd and 4th columns we give the values for Ar^ = r^ p — r*^ obtained from our new fits. The 
last column contains the values for Ar c /j for the SLy7 Skyrme interaction taken from Ref. Q|. 



Nuclei 


exp 
r ch 




ch 1 c h 1 ch 




KDE0 


KDE 


SLy7 


16Q 


2.730 


-0.041 


-0.039 


-0.017 


40 Ca 


3.490 


0.000 


0.011 


0.020 


48 Ca 


3.480 


-0.021 


-0.008 


-0.015 


56 Ni 


3.750 


-0.018 


0.000 


-0.008 


88 Sr 


4.219 


-0.002 


0.019 




90 Zr 


4.258 


-0.008 


0.013 




208pk 


5.500 


0.011 


0.041 


0.002 



TABLE VIII: Critical density p cr , rms radii of the valence neutron orbits r v and spin-orbit splitting 
(S-O). The experimental values ( and the theoretical error a) used in the fit to determine the Skyrme 
parameters are taken as follows: for the p cr we assume a value of 2.5po ip = 0.5po)> t ne values of 



r v were taken from Ref. 



44 



451 ] (a = 0.06 fm) and the spin-orbit in 5 "Ni were taken from Ref. 



(a = 0.2 MeV). In columns 3 — 6 we give the results obtained from our new fits. 





Expt. 


KDE0 


KDE 


Per/ Pa 


2.5 


2.5 


2.1 


r v {vld b/2 ) (fm) 


3.36 


3.42 


3.41 


r v )(fm) 


3.99 


4.05 


4.03 


e n (2 Pl/2 ) - e n (2p 3/2 ) (MeV) 


1.88 


1.84 


1.81 


e p (2pi /2 ) - e p (2p 3/2 ) (MeV) 


1.83 


1.64 


1.63 
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TABLE IX: Comparison of the breathing mode constrained energies (in MeV) obtained for the 
KDEO and KDE interactions with the experimental data. 



Nucleus 


Expt. 


KDEO 


KDE 


90 Zr 


17.81 


17.98 


17.91 


116 Sn 


15.90 


16.42 


16.36 


144 Sm 


15.25 


15.53 


15.47 


208 pk 


14.18 


13.64 


13.60 



TABLE X: Results for the neutron skin, r n — r p (in fm), for all the nuclei considered to obtain the 
KDEO and KDE interactions. 



Nuclei 








KDEO 




KDE 


16Q 


-0.031 




-0.025 


24 Q 


0.510 




0.510 


34 SI 


0.189 




0.192 


40 CA 


-0.051 




-0.046 


48 CA 


0.158 




0.159 


48 NI 


-0.282 




-0.274 


56 NI 


-0.056 




-0.052 


68 NI 


0.175 




0.174 


78 NI 


0.287 




0.285 


88 SR 


0.095 




0.096 


90 ZR 


0.064 




0.065 


ioo SN 


-0.081 




-0.078 


132 SN 


0.220 




0.217 


208p B 


0.160 




0.155 
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TABLE XI: Single-particle energies (in MeV) for 40 Ca nucleus. 



Orbits 


Expt. 




KDEO 


KDE 


Protons 


ls l/2 


-50±11 




-39.40 


-38.21 


1P3/2 






-26.95 


-26.42 


lPl/2 


-34±6 




-22.93 


-22.34 


14/2 






-14.49 


-14.51 


2si/ 2 


-10.9 




-9.48 


-9.66 


14/2 


-8.3 




-7.59 


-7.53 


1/7/2 


-1.4 




-2.38 


-2.76 


Neutrons 


l s l/2 






-47.77 


-46.13 


1P3/2 






-34.90 


-33.92 


lPl/2 






-30.78 


-29.73 


14/2 






-22.08 


-21.66 


2si/ 2 


-18.1 




-17.00 


-16.78 


14/2 


-15.6 




-14.97 


-14.48 


1/7/2 


-8.32 




-9.60 


-9.58 


2p 3 /2 


-6.2 




-4.98 


-5.15 
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TABLE XII: Single-particle energies (in MeV) for 208 Pb. 



Orbits 


Expt. 




KDEO 


KDE 


Protons 


109/2 


-15.43 




-17.85 


-17.34 


107/2 


-11.43 




-13.77 


-13.39 




-9.70 




-11.37 


-11.23 


1^11/2 


-9.37 




-9.87 


-9.68 


2^3/2 


-8.38 




-9.43 


-9.30 


3si/ 2 


-8.03 




-8.67 


-8.62 


1^9/2 


-3.77 




-4.00 


-3.99 


2/7/2 


-2.87 




-2.78 


-3.00 


1*13/2 


-2.16 




-1.62 


-1.72 


3P3/2 


-0.95 




0.60 


0.26 


2/5/2 


-0.47 




-0.19 


-0.42 


Neutrons 


1^9/2 


-10.85 




-12.39 


-12.24 


2/7/2 


-9.72 




-11.60 


-11.64 


1*13/2 


-9.01 




-9.33 


-9.20 


3^3/2 


-8.27 




-8.67 


-8.77 


2/5/2 


-7.95 




-8.59 


-8.64 


3P1/2 


-7.38 




-7.54 


-7.65 


2^9/2 


-3.94 




-2.86 


-3.06 


1*11/2 


-3.15 




-1.65 


-1.69 


lj'l5/2 


-2.53 




-0.41 


-0.43 


3d 5 /2 


-2.36 




-0.43 


-0.64 


4si/ 2 


-1.91 




0.08 


-0.08 


2^7/2 


-1.45 




0.38 


0.20 


3^3/2 


-1.42 




0.56 


0.40 
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